Spin-charge Separated Solitons in a Topological Band Insulator 
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In this paper we construct a simple, controllable, two dimensional model based on a topological 
band insulator. It has many attractive properties. (1) We obtain spin-charge separated solitons 
that are associated with ir fluxes. (2) It suggests an alternative way to classify Z2 topological 
band insulator without resorting to the sample boundary. (3) When the n fluxes are dynamical 
variables, as in a correlated insulator with emergent gauge fluxes, these solitons are propagating 
bosonic excitations and their condensation triggers a phase transition into a planar ferromagnet. 



PACS numbers: 

There has been recent interest in novel varieties of band 
insulators which differ in subtle but essential ways from 
ordinary band insulators. The best known example is the 
Chern insulator [1] which breaks time reversal symmetry. 
Tight biding models of such insulator has been studied 
long back by Hofstadter. More recently, a tight binding 
for Chern insulator which has no net magnetic flux was 
proposed by Haldane.[2] Recent progress has focused on 
time reversal invariant insulators, where a natural gen- 
eralization of the Haldane model emerges on including 
spin orbit interactions [3]. These models have attracted 
considerable attention recently because of their relevance 
to the quantum spin Hall effect. [6] Their band structures 
are characterized by non-trivial Z 2 topological quantum 
number[3, 5], which differentiates them from ordinary 
band insulators, hence the terminology topological band 
insulators (TBIs) . Here we point out a remarkable prop- 
erty of these insulators - the presence of spin charge sep- 
arated solitons in the presence of ir flux - which allows 
for a bulk definition of the TBI. Spin charge separated 
excitations are of tremendous interest in condensed mat- 
ter physics. While spin-charge separation is common in 
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Here i labels the sites of the square lattice, \ld — 
{c\ a , cJjcr) is a two-component electron operator where a is 
the spin index and 1, 2 are the flavor (or 'orbital') indices. 
These play the same role as sublattice indices in the hon- 
eycomb model. Also, t, A are real hopping parameters 
while v is an onsite orbital splitting energy. This free- 
fermion Hamiltonian is invariant under the z-axis spin 
rotation (Rs z )- Moreover, since the hopping matrices 
of the up spin fermions are the hermitian conjugate of 



ID, in higher dimensions it is extremely rare, requiring 
the presence of novel quantum states with topological 
order [7]. The spin-charge separated solitons identified 
here are not propagating excitations since they are tied 
to external tt flux. However, if the flux itself is a dy- 
namical variable, then we show that these solitons are 
bosonic excitations and different type of solitons obey 
mutual semionic statistics. 

The model We will start from a topological band in- 
sulator (TBI) fermion model, and then couple it with a 
dynamical Z2 gauge field. TBIs are free fermion insu- 
lators. Their band structures are characterized by non- 
trivial topological quantum number [1, 3, 5], which dif- 
ferentiates them from ordinary band insulators. Two 
well known examples of TBI are the models proposed 
by Haldane [2] (time- reversal breaking) and the Z2 insu- 
lator model proposed by Kane and Mele (time reversal 
invariant) [3] on the honeycomb lattice. The Hamiltonian 
of our TBI, which is a square lattice version of the Kane- 
Mele model, is given by the Hamiltonian below with the 
bond variables rf set to +1 everywhere: 
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that of the down spin fermions, it is also invariant un- 
der time- reversal (T). So long as A ^ this model is 
in the TBI phase when < \v\ < 8\t\. Note that if one 
ignores the spin down fermions and replaces C2 by c{ for 
the up spin fermion, the free fermion Hamiltonian dis- 
cussed above becomes the Bogoliubov Hamiltonian of a 
spin polarized p+ip superconductor [8]. The main differ- 
ence between our free-fermion Hamiltonian and that of 
the p+ip superconductor is symmetry : while charge is 
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conserved in the former (hence possess global f/(l) sym- 
metry), it is only conserved modulo two (hence possess 
global Z 2 symmetry only) in the latter. It is well known 
that the vortex cores of the p + ip superconductor has 
Majorana fermion zero modes [8]. As shown in Ref.[9] 
the corresponding U(l) symmetric model describes a TBI 
and possess a pair of opposite spin, counter-propagating, 
edge states at the sample boundary, analogous to the 
Kane-Mele model [3]. 

In Ref. [10] it was shown that when the hopping be- 
tween the edges is reinserted with a twist in sign, equiv- 
alent to flip the sign of half a row of t^+j, in Eq. (1) 
and hence introducing a 7r flux, an edge Jackiw-Rebbi 
soliton is created[10]. Such a soliton is a point defect in 
two-dimensions, which possesses two fermionic (not Ma- 
jorana fermion) zero modes, one for each spin[10]. We 
call these soliton defects "fluxons". 

Spin-charge separated fluxons If in a plaquette 
Y\ijeO T ij = — 1' *his signifies a fluxon. We have per- 
formed numerical calculations and shown that for a wide 
range of the hopping parameters v, t, A the creation en- 
ergy of two fluxons are lower than the minimum energy 
for particle-hole excitations, i.e., the band gap. As shown 
in Ref. [10], there are two fermionic midgap modes local- 
ized on each fluxon, which are Kramers conjugate. Since 
the model (1) enjoys two independent particle hole sym- 
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where a =]' or a =j, 



these modes are precisely at zero energy (i.e. in the mid- 
dle of the gap). 

The occupation/unoccupation of these zero modes 
leads to an excess/deficit of 1 /2 fermion number per spin. 
The four different ways of occupying these zero modes 
(Fig.(l)) give rise to four different types of fluxons with 
the following quantum numbers: +±j. (charge 1, 

S x = 0); /_i t _|| (charge -1, S z = 0); / + | T (charge 
0, S z — |); f_i t + k\ (charge 0, S z — — |). The presence 
of these modes, and their quantum numbers, can also be 
deduced from flux threading arguments for the up and 
down spin integer quantum Hall states. In the absence 
of particle-hole symmetry the modes are no longer pre- 
cisely at zero energy, but must still be within the gap. As 
discussed subsequently, this structure is essentially pre- 
served even when spin rotation symmetry is completely 
broken, as long as time reversal symmetry remains. In 
Fig. 2 we present the charge and spin density profiles for 
a pair of + 1 j_ and _Aj fluxon. 

Quantum statistics of fluxons When fluxons are 
mobile their quantum statistics becomes important. We 
determine their statistics through explicit computation 
of the Berry's phase. The statistical angle between two 
fluxons (not necessarily identical), 0(/i,/2), is defined 
as 1/2 times the difference of the following two Berry's 
phases. The first is obtained by hopping /i in a clockwise 
loop enclosing / 2 , while the second is obtained by hop- 
ping /i along the same path but with / 2 sitting outside 
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FIG. 1: Four different types of fluxons. The notation, e.g., 
_i| means the up-spin zero- mode is filled while down- 
spin zero-mode is empty. 




FIG. 2: The charge (a) and spin (b) densities of a pair of 
fluxons on a 24 x 24 lattice with periodic boundary condition. 
The fluxon at coordinate (6,6) is It has charge 1 

and S z = 0. The fluxon at (18, 18) is it has charge 

and S z = 1/2. Note that that the average charge density 
of the ground state is 2. The spin-charge separation in real 
space is manifest. 



the loop. Given a closed loop sequence of fluxon posi- 
tions {(x^,x|), i = 0,l,2,---iV} the Berry's phase is 
given by 
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where H hop = JJ2i( T Lx + T i, y )i and l$(xi,x 2 )} is the 
fermion many-body ground state consistent with two 
fluxons being at Xi and x 2 [12]. Since the up-spin band 
and the down-spin band decouples, the whole electronic 
wavefunction is a product of two Slater-determinants 
|$(xi,x 2 )) = |$(xi,x 2 )) T ® |$(xi,x 2 ));. As a result, 

0(/aiTAl>/aaT,/93.1.) = ^(/aiT) /oat) + 0(fPil > ffal ) ' ( 3 ) 

where 0(/ ai f , fa 2 ]) IS th e statistical angle between the 
two fluxons in the up-spin band. The results for 
0(f ai f, /a 2 t) are presented in Fig. (3). They are con- 
sistent with the statistics obtained from anyon fusion 
arguments [11]: Let us discuss on the up-spin band only. 
Consider a bound state of two of fluxons /if, and an- 
other bound state of two fluxons /i | . Then each bound 
state carries charge 1 and flux 27r ~ and thus is a 
fermion. As a result the statistical phase between two 
/if fluxons would be one-quarter of that of fermions, i.e. 
±j. By numerical calculation we find 9(fip /i*) = — \. 
From particle-hole symmetry we immediately conclude 
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FIG. 3: The result for 0(/it,/2t) on a 4d x 4d square lat- 
tice with periodic boundary condition. We fix the posi- 
tion of /2f and let /if loops around along a 2d x 2d square 
path. The fermionic band parameters used in the compu- 
tation are: u — 0.3, A = 0.5, t = l.(a) Identical fluxons 
6 id = fl(/i|,/i f ) = 6{f_ i T , /_ and (b) distinguishable 
fluxons # op = 0(/_i|,/if) = Note that due to 

the particle-hole symmetry = if, /_ i ♦) and 

d(f_i+, /if) = #(/i |, The extrapolation to d — > oo 

gives -0 id = 6> op = 7r/4. 



that = —J, too. Now consider a bound 

state of an /i* fluxon and an fluxon. This bound 

state carries charge and should be a boson. This implies 
that = f. Oifcat, fct 2 t) one can determine 

the statistical phase in the down-spin band 9(f ai ±, f a2 i) 
readily: 
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This is because the Hamiltonian for the down spin band 
is the hermitian conjugate of that for the up spin band. 
Given Fig. (3) and Eqs.(3,4) we have determined the 
quantum statistics of fluxons. The result is shown in 
the following Table. In general fluxons should experi- 
ence a background magnetic field (the fermion density) 
as they hop around. However, since there are on aver- 
age two fermions per site (see Fig. 2(a)), this background 
magnetic flux is 27r per plaquette, hence is equivalent to 
no flux. The above results should be robust against per- 
turbations so long as the bulk gap is preserved. 



Self statistics 


Mutual statistics 


c = ±1/2 


c= ±1/2 


0(/cT,cj)/cT,cjJ = 


0(/cT,cj.> /-cT,-ci) = 


0(/cT,«4i /cT,-ci) = f 


0(/cT,c|.) /-cT,c|) = f 


0(/cT,-cj.j /cT.-cj.) = 


fl(/cT,-oi,/- C T,ci) = 


0(fcl,-cl, fchcl) = f 


0(/cT,-cl.) /-cT,-ci) = f 



A new way to diagnose Z2 TBI Note that the four 
fluxon states in Fig.l are degenerate due to T, PH^±. The 
degeneracy between the charged and neutral fluxon can 
be easily removed by adding a weak short range charge 
repulsion to the original fermion model. After that, one 
expects the lowest energy fluxons to be the neutral ones: 



/+|t -ix an< ^ ^-jT.+M - ^ n * ne res ^ °f ^ ue P a P er we refer 
to them as spin fluxons. The S z = ±1/2 spin fluxons 
form a Kramer's pair upon time reversal. 

So far in our discussion Rs z is a symmetry of the 
Hamiltonian. This global f (1) symmetry justifies the 
corresponding TBI to be called a U(l) TBI. How- 
ever, the presence of a Kramer pair of neutral fluxon 
is more general. We have checked that as long as 
T is unbroken, each neutral fluxon always comes as a 
Kramer pair. This is true even after breaking Rs z (by 
adding, say, T-invariant spin-flip hopping term to the 
TBI Hamiltonian [3]), and/or Piffj. (by adding, say, a 
chemical potential term to the TBI Hamiltonian). This 
robust degeneracy allows one to diagnose the T-invariant 
TBI, or Z2 TBI[3] without resorting to edge states. For 
example, consider a Z2 TBI on a torus. One can intro- 
duce 2N far apart, low-energy, spin fluxons by, e.g., im- 
posing an energy penalty for charge accumulation. The 
ground state will be 2 2Ar -fold degenerate. On the other 
hand a trivial band insulator has no such degeneracy. 
Hence this degeneracy differentiates a Z2 TBI from a 
trivial band insulator. This can be implemented as a 
numerical diagnosis of Z2 TBIs. 

This study naturally generalizes to three dimension. 
For the 3D-Z2 insulator (which is refered as the strong 
topological insulator in literatures, for instance [4, 5]), 
we find for a closed 7r-fhix loop, there are two gapless 
one-dimensional Dirac fermion modes propagating along 
the 7r-fiux loop in opposite directions and are Kramers 
conjugates of each other. 

Dynamical 7r-fluxes and TBI* In order to make 
the fluxon elementary excitations, we give the Z2 variable 
dynamics. This is achieved by adding the following term 
to Eq. (1). 
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The fermions in the above Hamiltonian carries a 
Z2 gauge charge, hence are not ordinary electrons. We 
refer to such a correlated band insulator with emergent 
Z2 gauge fields as a TBI*. Nonetheless, the fundamen- 
tal fermion degrees of freedom of Eq. (5) possesses both 
the fermion number and the spin quantum number. In 
the following we show that the elementary excitations of 
this model exhibit separation of the the fermion quantum 
number (which we abbreviated by "chanrge") and spin. 

In Eq. (5) the term n^eD 7 ^ is tne ^ 2 g au S e nux 
going through a plaquette, and K, J are gauge cou- 
plings. The last term of Eq.(5) causes the fluxons to 
hop from one to a neighboring plaquette. As usual, t x,z 
are the first and third components of the Pauli matri- 
ces. The Hamiltonian in Eq.(5) has to be supplemented 
with a local constraint on every site (the 'Gauss Law') 

Iljen.n.of i T i] = (-1)*'"* 1 '", where the product is over 
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nearest neighbors of the site 'i'. For J = the ground 
state of Eq.(l) lies in the gauge sector where there is no 
flux in any plaquette. Under that condition it is always 
possible to tune the parameters so that the fluxons are 
the lowest energy excitations in the fcrmionic sector. For 
non-zero J the static fluxons arc no longer eigen exci- 
tations. However, so long as the fluxon creation energy 
>> J the delocalization of fluxons will not close the ex- 
citation gap. In that limit the gapped mobile fluxons 
exhibit spin-charge separation as illustrated in Fig.l. 

Spin fluxon condensation In the rest of this paper 
we will assume the spin fluxons to be the lowest energy 
excitations. Now let us ask what happens as the magni- 
tude of J is increased. When the energy cost in creating 
a static spin fluxon is counter balanced by the kinetic 
energy gain due to its delocalization, spin fluxons will 
spontaneously proliferate. Owing to their Bose statistics 
this will trigger Bose condensation at zero temperature. 
It is interesting to ask what is the nature of the new 
ground state and what is the nature of the (quantum) 
phase transition. In the following we shall discuss two 
scenarios. 

(I) If S z is conserved, two spin fluxons of opposite S z 
can be created and annihilated dynamically, while two 
fluxons with the same S z can not. In this case we can 
view the S z = — \ fluxon as the anti-particle of S z = +\ 
fluxon, and T transforms one into the other. The sym- 
metry which dictates the S z conservation is Rs z - Un- 
der such condition, the field theory describing the spin 
fluxon condensation is characterized by the following La- 
grangian density 

£ = \\dA\ 2 + i|v^| 2 + + 1 U |0| 4 , (6) 

where <f> is the complex fluxon field. The two phases of 
this field theory are: 1) the fluxon uncondensed phase 
where (0) = and Rs z is unbroken. In this phase, creat- 
ing a spin fluxon costs a finite energy. In the gauge the- 
ory jargon the Z 2 gauge field is in the deconfined phase. 
This is the phase of a spin liquid with a finite gap for 
spinon (bosonic) excitations. 2) The fluxon condensed 
phase where (</>} ^ and Rs z is spontaneously broken. 
This is a phase where the Zi gauge field fluctuates so 
strongly that it confines the fcrmionic charge excitations. 
Magnetically it is an XY ordered ferromagnet. (We have 
implicitly assumed that the ordering is easy plane rather 
than easy axis, which is natural in the presence of spin- 
orbit coupling [15]) Moreover, since the fermionic charge 
excitation are absent at low energies throughout the tran- 
sition, this phase is an electric insulator. Thus spin fluxon 
condensation triggers a spin liquid to a ferromagnetic in- 
sulator transition. According to Eq. (6), the universality 
class of the transition is 3D XY. The fact that <fi trans- 
forms as e 10 / 2 while the order parameter S + transform 
like e td under R$ z implies the identification S + ~ <p 2 . 
Hence, there is a subtle difference from the regular XY 



transition obtained from magnon condensation (i.e. con- 
densing S + itself), in that the order parameter's critical 
scaling dimension is anomalously large [13]. 

(II) S z is not conserved, but T is preserved. Now, 
one can add spin rotation breaking terms to the effective 
Lagrangian as long as they preserve time reversal sym- 
metry. The first such term in the long wavelength limit 
is g(S + ) 2 + h.c, is actually a quartic term when written 
in terms of the spinon fields introduced above gcj) 4 + h.c. 
Now, the condensation of <f> leads to a confined insulator 
with the spontaneous breaking of time reversal symme- 
try. Interestingly, although such an insulator has an Ising 
order parameter, the transition is expected to remain 3D 
XY like, due to the irrelevance of four fold anisotropy at 
the XY critical point. 

In the past, the transition to magnetically ordered 
states from spin liquids has been described using the 
Higgs mechanism. Here, we have described how confine- 
ment can also lead to magnetic order. This mechanism 
can lead to novel quantum phase transitions complement- 
ing those discussed in [14], which will be described in 
future work [15]. 

After completing this work, we learnt that in a re- 
cent preprint arXiv:08010252 X-L Qi and S-C Zhang have 
obtained similar results[16]. We thank Joel Moore and 
Cenke Xu for helpful discussions. The authors were sup- 
ported by the Directior, Office of Science, Office of Basic 
Energy Sciences, Materials Sciences and Engineering Di- 
vision, of the U.S. Department of Energy under Contract 
No. DE-AC02-05CH11231. 
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